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We apply the method of perturbative fragmentation to study bottom fragmentation in top 
decay. We present the energy spectrum of fo-quarks and b-navoured hadrons in top decay. 
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A reliable understanding of bottom-quark fragmentation in top-quark decay (t — > bW) will 
be fundamental to accurately measure the top properties, such as its mass, at present and future 
high-energy colliders. At the LHC, final states with leptons and J/ip, where the leptons are W- 
decay products and the J/ip comes from the decay of a ^-flavoured hadron, will be a promising 
channel to reconstruct the top mass nit ™- The error is estimated to be Amj ~ 1 GeV and 
the 6-fragmentation is the largest source of uncertainty, accounting for about 0.6 GeV. In this 
paper we investigate bottom fragmentation in top decay within the framework of perturbative 
fragmentations and use phenomenological models to describe the 6-quark hadronization. 

According to the factorization theorem, the rate for the production of a 6-hadron in top 
decay can be written, up to power corrections, as the following convolution: 
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m t , /jl f ) <g) D p (x b , m b , fx F ) <8> D np (x B ). 



(1) 



In Eq. (||), To is the Born width of the process t — ► bW, dT/dx b is the rate for the production 
of a massless b quark in top decay (coefficient function), D p (x b , m b , fj,p) is the perturbative 
fragmentation function which expresses the transition from a massless b into a massive b quark, 
Hf is the factorization scale, D np is the non-perturbative fragmentation function. x b and xb are 
the normalized 6-quark and 6-hadron energy fractions in top decay 
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We have computed the rate for the production of a massless b quark in top decay in dimen- 
sional regularization and subtracted the collinear singularity in the MS factorization scheme. 
We have got the following MS coefficient function: 
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The value of the perturbative fragmentation function D p (x b , m b , fj,p) at any scale \if can be 
obtained by solving the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equa- 
tions, once an initial condition at a scale [Iqf is given. The initial condition of the perturbative 
fragmentation function reads!: 



D p (x b , hq F , m b ) = 5(1 - x b ) + 
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The result in Eq. (Q) is process independent § and can hence be used in top decay as well. 
Moreover, it has been shown!, by comparing the rates for massless and massive 6-quark pro- 
duction in top decay, that the coefficient function (||) is indeed consistent with Eq. (|4|). Solving 
the DGLAP equations for the evolution //of — ► /J>f, with a NLO kernel, allows one to resum 
logarithms log(fi F / fJ^p) with NLL accuracy. If we set Hof — fn b and \xf — m t , we resum large 
logarithms ~ log(m|/m|) which appear in the massive, unevolved, fixed-order calculation of 
dT I dx b 0. The coefficient function (||) and the initial condition (||) contain terms which get large 
once x b — > 1. This corresponds to soft-gluon radiation in top decay. Soft terms in the initial 
condition of the perturbative fragmentation function are process independent and have been re- 
summed with next-to-leading logarithmic accuracy!. Soft-gluon contributions in the coefficient 
function are process dependent and their resummation is currently in progress. 

In Fig. |H we present the 6-quark energy distribution in top decay according to the approach 
just described. We note that the inclusion of soft-gluon resummation in the initial condition of 
the perturbative fragmentation function smoothens out the distribution and the x b spectrum 
exhibits the so-called Sudakov peak. Also shown is the massive 0(as) result, which lies below 
the other two distributions throughout the spectrum and is divergent as x b —*■ 1. It has also been 
found that the implementation of soft resummation yields a milder dependence of observables 
on the factorization scale hof which enters Eq. 

We would like to make predictions for the spectrum of 6-flavoured hadrons in top decay. For 
this purpose we use some phenomenological models which contain tunable parameters which are 
to be fitted to experimental data. We consider a power law with two parameters: 
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Figure 1: fa-quark spectrum in top decay according to the perturbative fragmentation method, with (solid) and 
without (dashes) NLL soft-gluon resummation in the initial condition (Q), and according to the O(as) massive 
result (dots). In the inset figure, the same distributions are plotted on a logarithmic scale. 



Table 1: Best-fit values for the parameters contained in the non perturbative fragmentation functions. 





ALEPH 


SLD 


a 


0.31 ±0.15 


1.88 ±0.42 


(3 


13.21 ± 1.62 


27.04 ± 4.02 


X 2 (a,/3)/dof 


2.62/14 


11.12/16 


<5 


20.39 ± 0.77 


18.80 ±0.60 


X 2 (5)/dof 


17.27/15 


17.46/17 


e 


(1.12 ±0.16) x 10~ 3 


(1.17 ±0.10) x 10~ 3 


x'(e)/dof 


22.96/15 


130.80/17 



the model of Kartvelishvili et al. 
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np (x;6) = (1 + 5)(2 + 5)(1 - x)x d 

and the non-perturbative fragmentation function of Peterson et al. 
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In Eq. (|5|), B(x,y) is the Euler Beta function; in (Q) A is a normalization constant. We tune 
such models to e + e~ data from the SLD I and ALEPrli Collaborations. The SLD data refers 
to 6- flavoured mesons and baryons and the ALEPH to mesons. When we do the fits, we must 
describe the e + e~ — > bb process within the same framework as for top decay, i.e. we use the 
perturbative fragmentation method, NLL DGLAP evolution and NLL soft-gluon resummation 
in the initial condition (Q). 

In Table || we report on the parameters which yield the best fits of the hadronization models, 
along with the corresponding standard deviations and the \ 2 P er degree of freedom. We note 
that the model in Eq. (||) fits best the data, although the best-fit parameters show pretty-large 
uncertainties. The Kartvelishvili model well reproduces both ALEPH and SLD data, while the 
Peterson model marginally agrees with the ALEPH data and is inconsistent with the SLD one. 
Moreover, the values of the parameters 8 and e, fitted to ALEPH and SLD, are in agreement 
within two standard deviations. In Fig. |2| we show our prediction for the 6-hadron spectrum 
in top decay, using all three hadronization models, fitted to ALEPH. We note that each model 
yields a statistically-different prediction for the xb spectrum, the Peterson-based distribution 



Figure 2: 6-hadron spectrum in top decay, describing the hadronization according to the power law with two 
parameters (solid), the Kartvelishvili (dashes) and the Peterson (dots) models, fitted to the ALEPH data. For 
any model, we plot the edges of a band corresponding to one-standard-deviation confidence level. 




Figure 3: As in Fig. g, using the power law, but fitted to ALEPH (solid) and SLD (dashes). 

being peaked at larger xb values. In Fig. || we compare the predictions obtained using the 
power law with two parameters, but fitted to ALEPH and SLD. We observe that the spectra 
are distinguishable; such a difference may be related to the different hadron types that the two 
experiments have reconstructed. 

In summary, we have reviewed the main results obtained once we studied the fragmentation 
of b quark in top decay using the perturbative fragmentation function approach and presented 
the energy distribution of b quarks and ^-flavoured hadrons in top decay. 
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